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Abstract 

We present an approach for self-consistent calculations of the many-body Green function in tran- 
sition metals. The distinguishing feature of our approach is the use of the one-site approximation 
and the self-consistent quasiparticle wave function basis set, obtained from the solution of the 
Schrodinger equation with a nonlocal potential. We analyze several sets of skeleton diagrams as 
generating functionals for the Green function self-energy, including GW and fluctuating exchange 
sets. Their relative contribution to the electronic structure in 3d-metals was identified. Calcula- 
tions for Fe and Ni revealed stronger energy dependence of the effective interaction and self-energy 
of the d-electrons near the Fermi level compared to s and p electron states. Reasonable agreement 
with experimental results is obtained. 
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Density-functional theory (DFT), in particular the local-density approximation (LDA), 
has proven to be a rather successful ab-initio approach to describe physical properties of 
many materials. Nevertheless numerous applications of this method have revealed a num- 
ber of shortcomings related to the inadequate treatment of both excited (energy gap) and 
strongly correlated states. For a long time the many-body Green function (GF) approach was 
considered as a possible ab-initio alternative to DFT, but its current applicability is usually 
restricted to the homogeneous electron gas (HEG) model fl], ^ || or semiconductors in the 
GW approximation. However, realization of the full self-consistent GW scheme is compli- 
cated due to enormous computational difficulties. As a result, for the transition metals (TM) 
the GW approach was applied only to Ni@, [5| and these calculations used the wave functions 
obtained within the LDA as the basis set, i.e. they were not self-consistent. Thus, it is quite 
unlikely that a universal GF method can be formulated and tailoring the approximations to 
specific properties of a given class of materials is desirable. The localized character of the 
d-wave function in TM often implies the well-known one-site approximation which in real 
systems acquires certain material specific features. 

Already the calculations of Refs. 0, |] revealed a strong energy-dependence of the effective 
interaction in TM, even in the vicinity of the Fermi level. For example in Ni the effective 
Coulomb interaction W is nearly zero at 6 eV and increases to an unscreened value of 20- 
25 eV at 20 eV. Another important observation is related to a proper intra-atomic scale 
treatment. Elementary analysis shows that the d- function in bcc Fe is significantly altered 
at a distance as small as 0.8 a.u., while s and p wave functions are altered at distances of 
about the Wigner-Seitz radius (~ 2.6 a.u.). The "local" screening length, A = xT 1 (>c 2 (r) = 
47re 2 n (p(r')), with n (p(r)) = (3/np) 1 ^ 3 /n), which varies from 0.5 a.u. in the region of 
the maximum of the d— wave function to a value of 1.0 a.u. at the atomic sphere (AS) 
boundary, provides an initial justification for the one-site approximation and suggests that 
the correct material-specific implementation must include a spatially resolved representation 
of the polarization operator (PO) and other two-channel operators, at least inside the AS. 
We believe that the energy dependence and intra-atomic resolution are both crucial for the 
quantitative description of TM and must be included explicitly in any reliable technique. 
On the other hand, we will show below that this approximation can describe the essential 
physics of TM, greatly reducing the computational efforts. 

In the present paper we incorporate the effects described above into the self-consistent 
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GF technique using the quasiparticle wave function basis set and the Luttinger-Ward 
functional]?! (LWF) approach for the self-energy calculations. We tested several sets of 
diagrams, mainly from the fluctuating-exchange subset (FLEX)||. All diagrams beyond 
these sets were taken in the local approximation. The use of the LWF guarantees the equiv- 
alence of one-particle properties calculated with the GF and with the corresponding total 
energy variation|J. The LWF formalism naturally leads to the integration over the imagi- 
nary axis, greatly improving the numerical accuracy of the integration. On the imaginary 
axis both the GF and PO quickly reach their quasiclassical limits, are sufficiently smooth, 
and are determined by the local potential only ||10|| . The rotation of the integration contour 
leads to a very convenient separation of the structural and local density dependent degrees 



of freedom for both self-energy and total energy |Tl|, pq| . Also, in contrast with the numerical 
method|| we analytically select the contributions from both the GF cut on the real axis 
and the integration along the imaginary axis. An energy linearization of the quasiparticle 
wave function similar to the one in linear density functional methods was used. 
In terms of an exact GF the thermodynamic potential is written as|7j 

n = -Tr {In [E - G l ] } - TrEG - $ (1) 

where E is the self-energy, G and Go are exact and 'bare' GF and $ is the Luttinger 
generating functional, which is represented by the set of skeleton graphs. Minimizing Q over 
G one can obtain 

The expressions for $ for the most important two-particle (hole) channels are summarized in 
Ref.|| The set from Fig. lb with 'bare' loops corresponds to the GW-approximation. After 
variation of $ over the GF, which we consider as a variational variable, one can obtain the 
usual GW expression for the self-energy: 

E( £ ) = - / G(e - u)VJl{u)W(u>)p: (3) 

where V c is the Coulomb interaction and W is the effective interaction W = V c /(1 + V C U). 
It is convenient to rotate the contour of integration in Eq.(||) in the complex plane, using 
the fact that H(z) has no singularities in the first and third quadrants[|l3|, [14]] , whereas G 
has a singularity (cut) in the third quadrant for e < Ep and in the first one for e > Ep. In 
the case of the quasiparticle the singularities of the GF are simple poles. After rotation the 



only additional contribution comes from the cut of the GF (Fig le), and for e < Ep this 
contribution is 

YP ca = J g(u')U{e-uj')V c W{£-uj')duj' e<0, 

where g(u) = — ImG(u;)sgnu;/7r. For the term corresponding to the integration along the 
imaginary axis we obtain 

oo oo 
-oo 

At Ep, e = and the pole term disappears. The term S* a is usually negative in the HEG 
model. On the real axis II (a>) has an imaginary part, which is important for accurate calcu- 
lations of W{oS) in TM. The summation over k (one-site approximation) transforms all the 
above formulas into matrix equations on the local wave function basis, greatly reducing the 
computational efforts. The full on-site GF G (r, r', e) = J2k (Gq 1 (r, r', e, k) — £ (r, r', e)) 
was obtained self-consistently from Eq.(|3|). Due to the large value of dTj/de in TM, the final 
G differs significantly from the initial Gq. The quantity II (r, r', e) = f G (e) G (e + uj) de/2m 
was also obtained self-consistently with the full G. 

Let us discuss the choice of skeleton graphs for the LWF. According to perturbation 
theory with Coulomb parameter a = e 2 m/pf the exchange diagram (Fig la) gives the largest 
contribution. The set of empty bubble diagrams (Fig. lb) gives the next term a 2 In 1/a. The 
summation of these terms is necessary to take into account the long-distance character of 
the Coulomb interaction. In addition, the proper treatment of the "dressed" GF in the 
bubble approximation must include vertex corrections |Tj| (Fig. lb). The next term (Fig.lc) 
is just the exchange diagram in the second-order approximation, which has an order of a 2 . 
It can be important in magnets, as the corresponding self-energy depends only on the GF 
with the same spin, while in bubble diagrams both spins are averaged. In principle, the 
Fig.lc diagram is already included in the set, if the bubble diagrams with vertex corrections 
are considered. But it seems desirable to include also the ladder sequence of Fig. Id, as it is 
also highly spin dependent. This set does not contain the first (a 2 ) term||. Such a sum is 
reduced by the effective screening of the Coulomb potential, which in turn strongly depends 
on the energy. The static value of the effective interaction is rather small, but it quickly 
increases as a function of energy (Fig.2b). We evaluated this sum of lattice-type T-matrix 
diagrams (Fig. Id) with the effective interaction replacing the energy dependent potential 
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by the averaged static interaction (V c . The parameter ( was chosen in such a way that 
the value of the diagrams on Fig.lc obtained with both energy dependent and (V c type of 
interactions were the same. The value of ( turned out to be ~ 0.35. 

The diagrams discussed above are also the leading diagrams from the point of view 
of the one-site approximation. The number of d-electrons with various angular moment 
projections m z is approximately the same in both Ni and Fe. If we take into account 
that the main contribution in the Coulomb interaction, expanded in spherical harmonics 
V(r — r') = ~^2 L Yj u (r)Yu(r t )(r < /r > ) is from the term with I = 0, then, with this accuracy, 
m z is conserved and one can classify various diagrams with the parameter 1/N = 1/(21 + 1). 
For d electrons this is a small parameter and the main terms with 1/N are again FLEX 



diagrams flBl- In summary, we believe, that the set of diagrams in Fig.l is the minimal set 
which must be included in the calculations. On the other hand, our calculations revealed 
that it can be sufficient for description of magnetic properties and electronic structure of 3d 
metals. 

The rotation of the integration contour allows us to deal only with the states lying in the 



vicinity of -£^[12], as the PO is smooth and decreases when Im e is increased, contrary to its 
behavior on the real axis. This is shown in Fig.2a,b for ferromagnetic Fe. For the states near 
Ep we can use the GF constructed with low lying LMTO states \I/^(r) = J2l a k a 0r( r )^i(^)> 
where $(r) and (fi](r) are the corresponding solutions of equations 

95V?(r) = (s - f )$(r) - f V(r, r')$(r')d 3 r' = (4) 



G^Kr) = -tf (r) - / 9E{£ ^ r 'h Kr')d 3 r', (5) 

with V(r, r') = E(ej, r, r'). E\ is the center of gravity of the band with orbital moment / and 
the eigenvectors of the generalized eigenvalue problem with the Hamiltonian 

H = H + - S,(eO - (e - efitfa) (6) 

which is a matrix in the space of LMTO states <fi + h<f), where Hll' are the usual LMTO 
coefficients. This procedure takes into account the energy-independent exchange potential 
H x = l/2m J V c lmGde exactly. 

The important problem in the nonlocal calculations is the inclusion of the valence-core 
interaction. In our method we included the core-valence exchange term exactly, whereas 
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the core- valence correlation term was added using the approximation developed in Ref.[17 



The latter term in TM such as Fe, Ni or Cu is small. The set of Eqs.(|]),([5]) was solved 
iteratively[|l^]. The integral equation for the screened interaction D(u,r,r') was solved 
using the product basis introduced in Ref . JT^j . 

The procedure above was applied to the 3d transition metals Fe and Ni. The main 
contribution to the self-energy comes from the exchange diagram (Fig la), but the pure 
exchange approximation produces too large a self-consistent magnetic moment M (M ~ 



3.05/% in Fe (see also Ref. [flQfl )). The bubble diagrams screen the Coulomb interaction, 
either too strongly (M ~ 1.95/% in Fe) or too weakly (M ~ 0.73/% in Ni). Vertex corrections 
only slightly modify this result, uniformly reducing the interaction for both spins. The 
conclusions in general agree with GW results for the HEG[P, 0. H], though the importance of 



the exchange T- matrix diagram (Fig Id) is rather intrinsic to TM[|, 20 1. With all considered 
diagrams the calculated equilibrium magnetic moments (2.04/% in Fe and 0.65/% in Ni ) 
are close to LDA values (2.15/%, 0.62/%) and experimental values (2.08/%, 0.59/%). 

The GF calculations naturally provide valuable information about the renormalization 
and damping of the electron spectra which are absent in LDA. In Fig. 3 we present the 
energy dependence of the real (S R ) and imaginary (T, 1 ) parts of the self-energy E 2a - ^ R (u;) 
is linear in a wide range of energies, supporting the concept of well-defined quasiparticles. 
On the other hand the renormalization factor Z\ = 1/(1 — dT^/de) strongly depends on 
the orbital number I. For s and p electrons Z\ is about 0.96 and is in good agreement 
with the HEG estimations {Z ~ 1 — 0.04e 2 m/p/ fl3||). But for d electrons this factor is 
about 0.6 in Ni and 0.7 in Fe, pointing out severe many-body effects. The value Zd was 
calculated self-consistently and in Ni it approximately coincides with the value calculated 
with a non-self-consistent (using LDA wave functions) method^. 

Previous studies revealed that the most significant differences between the LDA and 
experimental results are the satellite peak below the quasiparticle band (approximately 6-8 
eV below Ep) and a larger value of the density of states (DOS) at the Fermi level N(Ep). As 
a result, the overall experimental bandwidth (7-8 eV in Fe and 6-7 eV in Ni) is larger than 
what is found in the LDA along with the larger values of N(Ep). Extracted from the linear 
term in the heat capacity, the experimental DOS values are 57 st/a.u. in Fe and 81 st/a.u. 
in Ni. Though a considerable amount of this value can be ascribed to the electron-phonon 
interaction, they are still much larger than the 'bare' LDA values in Fe (30 st/a.u.) and in 
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Ni (48 st/a.u.). 

The small value of Zd obtained above leads to a narrowing of the quasiparticle band. The 
rest of the one-electron states take part in the formation of the satellite structure, which 
usually appears when the hole- hole interaction is taken into account pU| . Consequently our 
DOS (Fig. 4) is also wider compared to the LDA DOS and has a bump below the quasiparticle 
DOS. The value of the DOS at the Fermi level is even lower than in the LDA, though the 
value of coefficient 7, which is 1/Z times larger than the DOS, is approximately the same 
as in LDA. Both values are still lower than the experimental values. It can be due to the 
high sensitivity of the Stoner splitting in ferromagnets (and in turn the value of N(Ep)) to 
all the approximations used or possible narrowing of the hopping between neighboring sites 
in the spirit of the Gutzwiller approximation or the Dynamical Mean Field Theory, which 
is omitted in our approximation. 

In conclusion, we have proposed a new self-consistent version of the GF approach, which 
uses the quasiparticle wave functions as a basis set. We found that self-consistent GW ap- 
proach produces the leading contribution to the electronic structure and magnetic properties 
of TM, whereas the addition of fluctuating exchange diagrams overall slightly corrects this 
result improving comparison with experiment. While the self-consistent renormalization 
factors Z\ for s and p electrons in Ni and Fe are close to the estimations obtained from 
the HEG, Zj, is much smaller (0.6-0.7). The values of bandwidth and density of states at 
the Fermi level are in reasonable agreement with experiment. The proposed technique can 
be naturally used for total energy calculations. In summary we believe that the proposed 
technique can be considered as a practical ab-initio alternative to modern DFT methods 
with much wider range of applicability. 
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FIGURE CAPTIONS 

Fig 1 The diagram representation of: a) exchange energy, (b) "bubble" set with vertex 
corrections for the polarization operator, (c) exchange diagram of second order (d) exchange 
T-matrix set. The integration contour and pole positions in the complex energy plane are 
shown in (e). 

Fig 2. The energy dependence of the real and imaginary parts of the polarization operator 
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with / = 2 on the real axis (a). The same for the polarization operator and effective 
interaction V e ff = DiZi/Z on the imaginary axis (b). 

Fig 3. The energy dependence of the partial components of self-energy in ferromagnetic 
Fe: (a) the real part and (b) absolute value of imaginary part. 

Fig 4. The density of states in bcc Fe for the majority (solid) and minority (dashed) spin 
states. 
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